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Abstract 

We compare the diffuse Ha map of our Galaxy with the FUV (1370-1710A) continuum map. The Ha 
intensity correlates well with the FUV intensity. The Ha/FUV intensity ratio increases in general with the 
Ha intensity and the FUV hardness ratio (1370-1520A to 1560-1710A), implying that late OB stars may 
be the main source of the Ha recombination line at high latitudes. The variation of the Ha intensity as a 
function of the Galactic latitude is also very similar to that of the FUV intensity. The results likely suggest 
that not only the original radiation sources of the Ha and FUV backgrounds but also the radiative transfer 
mechanisms responsible for the diffuse backgrounds are largely common. Therefore, we propose a scenario 
wherein the Ha background at high latitudes is mostly composed of two components, Ha photons produced 
by in-situ recombination at the ionized regions around late OB stars and dust-scattered light of the Ha photons 
originating from late OB stars. 

Subject headings: diffuse radiation — ISM: structure — ultraviolet: ISM 



1. INTRODUCTION 

The diffuse Ha background is generally believed to 
originate from the extended, nearly fully ionized regions 
(called warm ionized medium, or WIM) along the lines 
of sight (Reynolds et al. 1998; see Haffneretal. 2009 for 
a recent review). The volume filling fraction of the 
WIM increases from ~ 0.1 at the midplane to > 0.3-0.4 
at \z\ = 1 kpc (Kulkarni & Heiles 1987; Reynolds 1991; 
Peterson & Webber 2002). The large energy requirement 
of the WIM strongly suggests that O stars in the Galac- 
tic plane are the primary ionization source (Reynolds 1984, 
1990a). There have been numerous investigations about how 
the Lyman continuum (Lyc) photons can escape the imme- 
diate vicinity of the O stars (or the bright H II regions) in 
the Galactic plane and ionize the diffuse interstellar medium 
(ISM) (Miller & Cox 1993; Dove & Shull 1994; Seon 2009a; 
Wood et al. 2010). However, the question of how the Lyc pho- 
tons can propagate almost freely up to the scaleheight (~ 1 
kpc) of the WIM has not been fully resolved because the pho- 
toionization models require the presence of very low density 
paths through the ISM. 

In this paper, we explore two alternatives to explain the dif- 
fuse Ha emission of our Galaxy. The first scenario, which was 
investigated for external galaxies, is that late, field OB stars 
outside the bright H II regions ionize some of the WIM, reduc- 
ing the required Lyc leakage from the bright H II regions in 
the galactic plane (Hoopes & Walterbos 2000; Hoopes et al. 
2001). The mean spectral type of stars in the field is later than 
that in H II regions (Massey et al. 1995; Hoopes & Walterbos 
2000). As the Ha intensity is directly proportional to Lyc in- 
tensity in the ionization bounded regions and the later-type 
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stars would produce a lower ratio of ionizing Lyc to non- 
ionizing far-ultraviolet (FUV) luminosities, the Ha/FUV in- 
tensity ratio can constrain the spectral type of the ionizing 
stars. Hoopes & Walterbos (2000) and Hoopes et al. (2001) 
measured Ha/FUV intensity ratios for 10 nearby spiral galax- 
ies, and concluded that late OB stars in the field are indeed an 
important source of ionization in those galaxies. 

Secondly, we note that the majority of the Ha excess in- 
tensity in a number of high-latitude clouds in our Galaxy 
can be understood as light scattered off the interstellar dust 
(Mattila et al. 2007; Lehtinen et al. 2010; Witt et al. 2010). 
The possible presence of scattered Ha radiation at high- 
latitude clouds was first considered by Jura (1979). Pho- 
toionization by an external field of Lyc is expected to re- 
sult in ionization of a thin outer shell of an isolated inter- 
stellar cloud while the observations of the Ha light from the 
high-latitude clouds clearly showed the dense inner core of 
the cloud as the brightest feature. Witt et al. (2010) argued 
that both a substantial fraction of the diffuse high-latitude Ha 
background intensity and much of the variance in the high- 
latitude Ha background are due to scattering by interstellar 
dust of Ha photons originating elsewhere in the Galaxy. If 
this is the case, the diffuse Ha emission should be closely 
related to the far-ultraviolet (FUV) continuum background, 
which is mostly dust-scattered light of the FUV stellar radia- 
tion (Bowyer 1991; Seon et al. 2011). 

Therefore, comparison of the diffuse Ha emission with the 
FUV continuum background could not only provide a clue to 
the spectral type of ionizing sources of the WIM in our Galaxy 
but also confirm the significance of dust-scattering in the dif- 
fuse Ha. In fact, Seon et al. (2011) found a strong correla- 
tion between the diffuse FUV (1370-1710A) and Ha back- 
grounds, suggesting a similarity between the origins of both 
backgrounds. In this paper, we present a more detailed com- 
parison of the Galactic Ha map with the FUV continuum map 
and propose a scenario on the origin of the diffuse Ha emis- 
sion. We describe the data used in the present study in Section 
2. In Section 3, the spectral type of the stars responsible for 
the diffuse Ha emission is investigated. Section 4 contains a 
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discussion on the dust-scattering of Ha photons. We discuss 
the main result in Section 5. A brief summary is presented in 
Section 6. 

2. DATA 

The FUV data was obtained by the Spectroscopy of 
Plasma Evolution from Astrophysical Radiation (SPEAR) in- 
struments, also known as Far Ultraviolet Imaging Spectro- 
graph (FIMS) (Edelstein et al. 2006a,b; Seonetal. 2011). 
The SPEAR/FIMS instrument is a dual-channel FUV imag- 
ing spectrograph ("Short" wavelength channel [S-band] 900- 
1150A, "Long" wavelength channel [L-band] 1350-1750A; 
A/AA - 550) with a large field of view (S-band, 4° x 4'.6; 
L-band, 7°. 4 x 4'. 3) and 10' imaging resolution onboard the 
first Korean astronomical satellite STSAT-1, designed to ob- 
serve diffuse FUV emission. The SPEAR/FIMS survey obser- 
vations were performed by scanning the sky at constant eclip- 
tic longitude from the north ecliptic pole to the south ecliptic 
pole, during the eclipses. The pole-to-pole scanning drifted 
360° along the ecliptic equator for one year because of the 
properties of a sun-synchronous orbit. The mission observed 
~ 80% of the sky. 

The diffuse FUV map in L-band was obtained by eliminat- 
ing stars from the total FUV sky map and by binning pho- 
ton events using the Hierarchical Equal Area isoLatitude Pix- 
elization (HEALPix) tessellation scheme (Gorski et al. 2005) 
with ~ 1° pixels (resolution parameter V s id e = 64). We 
excluded the wavelength regions of the strongest emission 
lines (Si IV A1398, Si II* A1532, C IV AA1548, 1551, and 
Al II A 1671) and averaged the data in the wavelength re- 
gions of 1370-1520A, 1560-1660A, and 1680-1710A. A de- 
tailed description of the FUV continuum map is presented in 
Seon etal. (2011). 

The all-sky survey maps in Ha and 100 /iin were obtained 
from Finkbeiner (2003) and Schlegel et al. (1998), respec- 
tively. The all-sky map of H I column density is from the 
Leiden/Argentine/Bonn (LAB) Survey (Kalberla et al. 2005). 
The data interpolated onto a HEALPix projection were ob- 
tained from the Legacy Archive for Microwave Background 
Data Analysis (LAMBDA). 

3. SPECTRAL TYPE OF IONIZING SOURCES 

Figure 1 shows the Ha/FUV intensity ratio maps obtained 
by dividing the Ha intensity with the FUV continuum inten- 
sity and smoothed with a spherical Gaussian function with 
a full-width at half maximum (FWHM) - 3° (Seon 2006). 
The smoothing was performed using the spherical harmonic 
transform in the HEALPix scheme so to smear the same sur- 
face areas in all directions. We excluded the regions with 
low signal-to-noise ratio (S/N < 3) and two regions between 
edges of observing strips centered at (l,b) ~ (45°, -35°) and 
~ (225°, 45°), which appeared to produce an artifact with 
very high Ha/FUV intensity ratios after the smoothing. The 
FUV and Ha intensities are expressed in units of rayleighs 
(1 rayleigh [R] = 10 6 /47r photons cm" 2 s -1 sr" 1 ) and contin- 
uum unit (1 CU = photons cm" 2 s" 1 sr -1 A" 1 ). In the fig- 
ure, the region with higher Ha/FUV ratio is represented in 
blue. Top and bottom panels show the intensity ratio maps be- 
fore and after dust-extinction correction, respectively. Using 
the standard approach of Bennett et al. (2003) for the emit- 
ting medium coexisting with dust and the extinction map of 
Schlegel et al. (1998), the correction for dust-extinction was 
performed for both radiations. In the top panel of the figure, 



it is clear that bright H II regions have higher Ha/FUV ratios 
than the diffuse Ha regions. Correction for dust-extinction 
lowers the ratio, because the correction for FUV is larger than 
for Ha. As dust-extinction is stronger at low latitudes, the 
correction for H II regions is bigger than for the WIM. There- 
fore, the distinction between the bright H II and the WIM re- 
gions becomes less clear after the dust-extinction correction, 
as shown in the bottom panel. 

For the quantitative comparison, the WIM regions (we will 
use the term WIM to denote the diffuse Ha regions that are 
not related to the bright H II regions, regardless of physical 
origin) were isolated in the Ha image using a masking pro- 
cedure similar to that described in Hoopes et al. (2001). The 
Ha image was smoothed with a median filter using a pixel 
size of ~ 15° which corresponds approximately to the size of 
the largest H II regions. The smoothed image was subtracted 
from the original, leaving an image of the small-scale struc- 
tures only. In other words, the smoothly varying component 
was subtracted. On this image a mask was created by as- 
signing the value ' 1' to the pixels with intensity lower than 
a prescribed threshold and '0' to the others. We tried sev- 
eral threshold values between 0. 1 and 50 R as the separation 
between the WIM and H II regions was rather fuzzy. The un- 
smoothed Ha image was then multiplied by the mask, leaving 
an image of only the WIM. A WIM FUV image was also cre- 
ated using the same mask. The Ha and FUV images for H II 
regions were created by multiplying the inverted WIM mask. 

Figure 2(a) shows the histograms of the Ha/FUV intensity 
ratios for the WIM and H II regions in our Galaxy, with the 
thresholds of 1, 5, 10 and 20 R. The threshold '1 R' masks 
most of the H II regions and '20 R' the brightest ones. The 
mean ratios for the H II regions are generally higher than the 
mean ratios in the WIM. The difference in the logarithm of the 
average ratio between the H II regions and the WIM increases 
smoothly from 0.40 ±0.40 to 0.62 ±0.39 dex as the thresh- 
old increases from 1 to 20 R. These values are smaller than 
those found in 10 external galaxies in Hoopes et al. (2001), 
although the largest value of ~ 0.62 is consistent with the re- 
sults for NGC 2903, NGC 5457, and NGC 1512 within the er- 
ror ranges. Figure 2(b) shows the histograms of the Ha/FUV 
intensity ratios after dust-extinction correction. Correction for 
dust-extinction lowers the difference in the average ratio, as in 
the bottom panel of Figure 1 . 

The smaller difference of the mean ratios in our Galaxy 
than those of other galaxies may be attributable to geomet- 
rical effects, as the WIM and H II regions overlap in our 
Galaxy while in the face-on external galaxies examined in 
Hoopes & Walterbos (2000) and Hoopes et al. (2001) two re- 
gions are better separated from each other. The field stars 
in the Galactic plane overlapping with H II regions would 
lower the Ha/FUV intensity ratio observed toward the bright 
H II regions. Some of the Ha and FUV photons origi- 
nating from H II regions are scattered into high latitudes 
(Wood & Reynolds 1999; Witt et al. 2010; Seon et al. 201 1). 
The dust-scattering of Ha and FUV photons from the H II re- 
gions to high latitudes would then increase the ratio for the 
WIM, further lowering the difference of the ratios. Further- 
more, all the H II regions are at the same distance in a face-on 
external galaxy while in the Milky Way the apparent size of 
H II regions is potentially dominated by their distances mak- 
ing the separation more difficult. Therefore, the difference 
between the H II regions and the WIM in our Galaxy may not 
be significantly different from other galaxies when the Galaxy 
is viewed face-on. 
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Along the top axis of Figure 2 are shown the stellar spectral 
types that would yield the corresponding Ha/FUV ratios. We 
derived a temperature and effective gravity for given spec- 
tral types using the calibration of Vaccaetal. (1996). We 
then calculated the Lyc and FUV photon luminosities L Lyc 
and Lfuv for spectral types of B2-03 stars in the wavebands 
A < 912A and 1370- 1710A, respectively, by interpolating on 
a grid of Kurucz models (Castelli & Kurucz 2003). The Ha 
luminosities Lhq equal 0.46LLyc, assuming case B recombi- 
nation. Since Vaccaetal. (1996) covers spectral types only 
earlier than or equal to BO. 5, we extrapolated the calibration 
down to B2. The temperatures estimated in this way were 
~ 5 - 19% higher than the values obtained from the older cal- 
ibration of Straizys & Kuriliene (1981). The extraporation to 
B2 would be accurate within ~ 19%. Figure 2 suggests that 
the difference between the Ha/FUV ratios in the WIM and 
those in the H II regions is due to the difference in main ion- 
izing sources. 

An interesting result in Figure 2 is that the median value of 
the observed ratios in the WIM corresponds to a value pre- 
dictable from ~ 09 - BO stars with effective temperatures of 
- 35,000 K. Sembach et al. (2000) modeled the WIM as a 
combination of overlapping H II regions and found that typi- 
cal values of the line ratios [S II] A6716/Ha ~ 0.2-0.4, [N II] 
A6583/Ha -0.3-0.6, [Ol] A6300/Ha ~ 0.01 -0.03, [OIII] 
A5007/Ha < 0.1, and He I A5876/Ha < 0.02 observed in the 
Galactic WIM can be well reproduced by an ionizing spec- 
trum of a star with r e ff — 35,000 K; Domg6rgen& Mathis 
(1994) also found that the WIM observations best agree with 
?eff ~ 38,000 K. The results are consistent with the Ha/FUV 
ratios, thereby strongly supporting that the Ha emission orig- 
inates from late OB-type stars. The optical line ratios will be 
discussed more in Section 5. We also note that the Ha/FUV 
ratio has a well-defined limiting value corresponding to B2 
stars, meaning that the latest stellar type that is capable of 
ionizing the WIM may be about B2. 

The mean Ha/FUV intensity ratio after and before the dust- 
extinction correction as a function of Ha intensity is shown 
together with the standard deviation of ratios from the mean 
value for a given Ha intensity inverval in Figure 3(a). Un- 
certainty of the ratio toward an individual sightline is much 
smaller than the deviation of ratios. Red and black sym- 
bols represent the intensity ratio before and after the dust- 
extinction correction, respectively. There is a general trend for 
the average ratio to increase with the Ha intensity, although 
it scatters a lot. Both the general trend and large scatter are 
qualitatively similar to those found in the external galaxies in 
Hoopes et al. (2001). The scatter in the ratio tends to increase 
with the Ha intensity. At low latitudes and highest Ha inten- 
sities, the Ha radiation may be dominated by recombination 
lines originated from H II regions around early O-type stars, 
while the FUV background is from dust-scattering (Seon et al. 
2011). Therefore, a much wider range of variability in the 
Ha intensity and the Ha/FUV ratio is found. On the other 
hand, at high latitudes and low values of Ha intensity, much 
of the structure in the faint Ha background may have the same 
origin as that of the FUV background and the variation of 
Ha/FUV decreases. In other words, the Ha contribution from 
dust-scattering becomes increasingly important at higher lati- 
tudes. A similar feature was found in Figure 21 of Seon et al. 
(2011), where a linear correlation of the Ha intensity with 
the FUV intensity was found only in low intensities. In addi- 
tion, at lowest Ha intensities, only early B-type stars would 



contribute to both backgrounds and the ratio would then con- 
verge to a single value. The converged value corresponds to a 
spectral type of B0 - B 1 . 

Interestingly, line ratios [S II]/Ha and [N II]/Ha, and the 
Ha/FUV ratio show a similar trend with decreasing Ha in- 
tensity, except the line ratios increase while the Ha/FUV ratio 
decreases. Mean values of the optical line ratios increase with 
decreasing Ha intensity. The ratios also show large scatter as 
the Ha/FUV ratio does. The mean Ha/FUV ratio decreases 
rapidly below Ha intensity of ~ 0.7- 1 R with decreasing Ha 
intensity in Figure 3(a). The mean line ratios increase most 
dramatically below the same Ha intensity (~ 0.7- 1 R) (Fig- 
ures 6, 12, and 15 in Madsen et al. 2006). The similar trend 
in both the optical line ratios and the Ha/FUV ratio may be 
understood as an increasing contribution of late OB stars to 
the WIM as the Ha intensity decreases. 

We also plotted the hardness ratio in FUV wavelengths 
versus the Ha/FUV ratio as the FUV hardness ratio might 
also provide an additional clue on the spectral type of the 
ionizing source of the WIM. The hardness ratio is defined 
as the ratio of average intensities at AA1370-1520A and 
AA1560- 1710A, excluding the bright C IV, Si II, and Al II 
emission lines, as in Seon et al. (201 1). We found an overall 
correlation between the FUV hardness ratio with the Ha/FUV 
intensity ratio, as shown in Figure 3(b), where the standard 
deviation of the FUV hardness ratios from its mean value for 
a given Ha/FUV bin is also indicated. The mean FUV hard- 
ness ratio shows rather rapid decrease below the Ha/FUV — 1 
(R/10 3 CU), corresponding to the approximate mean value of 
the Ha/FUV ratio of the WIM. Seon et al. (201 1) found that 
a map of the hardness ratio of 1370-1520A to 1560-1710A 
band intensity shows that the sky is divided into roughly two 
parts, which are determined by the longitudinal distribution of 
OB-type stars in the Galactic plane. Therefore, the correlation 
between the FUV hardness ratio and a large portion, if not all, 
of the Ha/FUV ratio supports that the diffuse Ha background 
is produced by relatively late-type stars rather than H II re- 
gions in the Galactic plane. 

Statistically significant correlation coefficients of 0.9 and 
0.8 for Figure 3(a) and (b), respectively, were obtained. We 
also note that most of the Ha/FUV intensity ratios are so con- 
centrated to ~ 1 (R/10 3 CU), as can be noted in Figure 2, that 
the two-dimensional histograms and/or contours as in Figure 
4 show strong peaks at the Ha/FUV intensity ratio of ~ 1 
(R/10 3 CU). Therefore, the average intensity ratio versus the 
abscissa in Figure 3 was the best way to represent clearly the 
correlations. We note that the Ha/FUV ratio map (Figure 1) 
does not match perfectly with the FUV hardness ratio map 
(Figure 14 of Seon et al. 201 1). However, the maps generally 
correlate well with each other, except a few regions such as 
near (l,b) ~ (30°, 30°). 

Hoopes & Walterbos (2000) and Hoopes et al. (2001) com- 
pared the observed Ha/FUV ratios with the theoretical calcu- 
lations assuming that the FUV photons are emitted from unre- 
solved point sources in the observed regions. However, FUV 
photons are well scattered by interstellar dust grains to large 
distances and the diffuse FUV emission is wide-spread over 
regions where no point sources are present (Bowyer 1991; 
Witt et al. 1997; Schiminovich et al. 2001; Seon et al. 201 1). 
The FUV intensity used in the analysis of face-on galaxies 
would be produced by in-situ point sources as well as dust- 
scattered light in the regions, whereas the FUV data in the 
present study is dominated by dust-scattered light originating 



4 



SEON ET AL. 



from late O- and early B-type stars (Henry 2002; Seon et al. 
2011). We also note that the spectral type of dominant radi- 
ation source of the FUV continuum background is consistent 
with that inferred from Ha/FUV intensity ratio. 

It would be interesting to estimate how a large portion of 
the ionizing power required for the Ha background can be 
provided by late OB stars. If the Ha emission originates from 
ionized gas in a plane-parallel geometry, the hydrogen recom- 
bination rate rc per cm 2 of galactic disk in the vicinity of the 
Sun is given by 

r c w 4 x I0% a (b) x sin s" 1 cm" 2 , 

where In a (b) is the Ha intensity at galactic latitude b 
(Reynolds 1984). Using an approximate Ha intensity distri- 
bution I\i a {b) ss 1 .Ocsc \b\ R, which was obtained from rather 
limited observations, Reynolds obtained rc « 4x 10 6 s" 1 
cm" 2 . However, using the result of Hill et al. (2008) obtained 
from the full WIM observations, the distribution of Ha at high 
latitudes is fitted on average by Iu a (b) ~ 0.625 esc \b\ R (see 
also Dong & Draine 201 1). This leads to a lower recombina- 
tion rate of rc ~ 2.5 x 10 6 s" 1 cm" 2 . The Lyc luminosities per 
cm 2 of the Galactic disk estimated for the solar neighborhood 
are (10-29) x 10 6 and (0.5-1) x 10 6 photons s" 1 cm" 2 for O 
and B stars, respectively (Reynolds 1984). The fraction of Lyc 
luminosity of 09 and 09.5 stars in total Lyc luminosity due 
to O stars is ~ 7- 10% (Terzian 1974; Torres-Peimbert et al. 
1974). The Lyc luminosity from the stars later than or equal to 
09 would then be (1.2-3.9) x 10 6 photons s" 1 cm" 2 . There- 
fore, late OB stars are able to account for at least one-half of 
the required ionizing power. 

In fact, Elmegreen (1975) considered the intensity of the 
Ha background at \b\ < 30° and concluded that isolated H II 
regions surrounding early B stars could account for up to one- 
half of the observed Ha background in some regions of the 
Galactic plane. He predicted that near the Galactic poles, B 
star H II regions would account for a mean Ha intensity of 
^0.1 R. This result is consistent with our estimation. 

4. DUST SCATTERING 

The late OB-type stars may not have enough power to pro- 
duce the whole Ha background if the Ha photons originate 
purely from ionized gas around late OB stars, as noted by 
Reynolds (1984) and in Section 3. However, if dust-scattering 
plays an important role, the dust-scattered Ha photons orig- 
inating from the sources other than O stars in bright H II re- 
gions would contribute a larger portion to the Ha background 
at high-latitudes and the severe requirement of Lyc leakage 
from H II regions in the midplane may be potentially alle- 
viated. In this section, we examine how the contribution of 
dust-scattering could be significant in the Ha background. 

Galactic observables may generally correlate with each 
other because of the general property of a plane-parallel ISM. 
The correlation between a pair of Galactic quantities would 
be enhanced when their sources and radiative transfer mech- 
anisms are analogous. If a large portion of the diffuse Ha in- 
tensity is due to dust-scattering, the Ha intensity would show 
a better correlation with FUV intensity than that with other 
quantities, such as neutral hydrogen column density N(H I). 
On the other hand, if the radiative transfer mechanism for the 
diffuse Ha emission is not related with the dust-scattering, the 
correlation between the Ha and FUV intensities would not be 
improved. 

Figure 4 shows the correlation plots of Ha intensity with 



100 fim dust emission and H I column density N(H I) in the 
same forms of two-dimensional histograms and contours as 
in Figure 21 of Seon et al. (201 1), in which the FUV intensity 
is compared with 100 /im intensity, N(H I), and Ha intensity. 
Correlation coefficients estimated in logarithmic scale are also 
shown in the figure. We find that the correlation of Ha inten- 
sity with FUV (shown in Figure 21(c) of Seon et al. 201 1) is 
a bit better, at least qualitatively, than with the others ISM 
tracers shown in Figure 4. The correlation relations with 100 
/iin intensity and N(H I) are a bit curved and flaring, whereas 
the correlation with the FUV background is quite straight. 
However, the correlation coefficient between the Ha and FUV 
backgrounds is not significantly higher than the other correla- 
tion coefficients. 

The fact that the correlation coefficients show no signifi- 
cant differences may be caused by large scatter in the corre- 
lation relations and the general similarity in the latitude de- 
pendence of the data sets (the H I column density, 100 fim, 
Ha, and FUV intensities). Although the Ha emissivity and 
H I column density have different dependences on the gas 
density, their latitude dependences have the same trend with 
inverse-sin | b | if they arise in plane-parallel geometries. The 
inverse-sin | b j dependence of the data sets dominates the cor- 
relations in Figure 4. To reduce the latitude dependence and 
scatter, we therefore plotted the average intensities multiplied 
by sin \b\ versus sin \b\ within each of the A sin \b\ = 0.01 lati- 
tude intervals for N(H I), 100 ^m, Ha, and FUV intensities 
in Figure 5. In the figure, the curves were paired and ar- 
bitrarily moved to compare the shapes. The top two curves 
compare the N(H I) and 100 fim emission. The middle and 
bottom pairs show the variations of Ha and FUV intensities 
after and before dust-extinction correction, respectively. The 
N(H I) and 100 /im emission are well described by the plane- 
parallel or inverse-sin \b law, whereas Ha and FUV intensi- 
ties deviate from the law. The deviation of Ha and FUV in- 
tensities from the inverse-sin \b\ law is mainly attributable to 
the effect of dust absorption and scattering. A radiative trans- 
fer model of diffuse Galactic light in plane-parallel dust layer 
shows the similar dependence of /sin|fo| on sin|fc| as in the 
bottom curves of Figure 5 (van de Hulst & de Jong 1969; see 
also Wood & Reynolds 1999 and Seon 2009b). We also note 
that the inverse-sin \b\ equation of the Ha intensity described 
in Section 3 is still a good approximation because of the major 
dependence of the Ha intensity on inverse-sin \b\, especially 
at high latitudes, even / sin \b\ is not a constant in Figure 5. 

We note from the bottom two curves before dust-extinction 
correction that the dependences of the Ha and FUV back- 
grounds on the Galactic latitude are very similar at high lat- 
itudes (sin | & | > 0.3) (See also Figure 13 of Haffneretal. 
1999 and Figure 9 of Seon et al. 2011), whereas at low lat- 
itudes (sin \b\ < 0.3) they show a difference. Stronger dust- 
extinction in FUV would suppress the FUV intensity near 
the Galactic plane, as shown in low latitudes. We therefore 
corrected the dust-extinction as denoted by the middle two 
curves (red and black solid lines for Ha and FUV, respec- 
tively). The dust-extinction correction for both backgrounds 
was performed using the standard approach of Bennett et al. 
(2003). The variations of both backgrounds are remarkably 
similar to each other. The FUV radiation can be emitted even 
from late B-type stars, whereas the Ha emission requires ex- 
citation from earlier stars. Therefore, the scaleheight of the 
FUV sources is a bit higher than that of the Ha sources. On 
the other hands, ionized H II regions must be more extended 
than point sources for the FUV photons, and the dust-scattered 
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Ha originating from H II regions should be more extended 
than those from point sources. Accordingly, the dependences 
of both backgrounds on the Galactic latitude are very similar 
even at low latitudes, after extinction corrections are applied. 

In Figure 6, we also examine correlations of the diffuse Ha 
intensity with N(H I), 100 /im emission, and the diffuse FUV 
background after the removal of inverse-sin \b\ dependence. 
In the correlation plot, we used only the latitude range of 
sin |£> | > 0.3 to avoid the strong dust-extinction effect. Red 
and black symbols represent the data points after and before 
the dust-extinction correction, respectively. The figure also 
shows correlation coefficients. It is obvious that the correla- 
tion with the FUV background is now significantly stronger 
than the correlations with others. In Seon et al. (2011), we 
compared the FUV continuum background with other ISM 
tracers and found that the FUV intensity correlates better with 
the Ha intensity than with others. The two-photon emission 
in the WIM is unlikely to be the cause of the strong correlation 
between the Ha and FUV backgrounds (Seon et al. 201 1). We 
are therefore convinced of the result that a significant portion 
of the diffuse Ha emission at high latitudes is dust-scattered 
light originated from elsewhere in the Galaxy. 

We now estimate the fraction of dust-scattered Ha emis- 
sion. Seon et al. (201 1) derived two equations (Equations (7) 
and (8)) relating the FUV intensity to 100 /im and Ha emis- 
sions at high latitudes. Since the constant terms in the re- 
lations are often regarded as extragalactic background, they 
should be identical. The difference between the constant 
terms is indeed negligible. Therefore, the relation between the 
100 /im and Ha emission are given by /ioo M m ~ 2.88 (Ihu/R) 
MJy/sr. Witt et al. (2010) provides an empirical relation, 
which relates the scattered Ha intensity to the 100 /im emis- 
sion at high latitudes. Using these relations, we obtain an 
average ratio of the scattered Ha intensity to total intensity, 
« 0.37/Ha- Therefore, in general, about 37% of the 
diffuse Ha emission at high latitudes would be due to dust- 
scattering. Contribution of the Lyc leakage from bright H II 
regions may be ~13% on average if late OB stars produce 
about half of the Ha emission. However, these numbers 
should be referenced with caution. The scattered portion can 
be much higher or lower depending on location, as shown in 
Witt et al. (2010). Since the FUV background traces the dust- 
scattered light, we expect that the FUV to Ha ratio would be 
more or less proportional to the scattered fraction of Ha back- 
ground. We note that the scattered fraction of Ha background 
shown in Figure 5 of Witt et al. (2010) is more or less anti- 
correlated with the Ha/FUV ratio map of Figure 1, verifying 
the expectation. 

5. DISCUSSION 

By combing the previous results, we propose that most of 
the diffuse Ha background at high latitudes may originate 
from relatively late O- and early B-stars outside the bright 
H II regions and from dust-scattered Ha photons from else- 
where in the Milky Way. The observed Ha/FUV intensity ra- 
tio in general increases with the Ha intensity, indicating that 
late OB stars could be important contributors to the diffuse 
Ha. As suggested by the strong correlation between the dif- 
fuse Ha and FUV backgrounds, dust-scattering may widely 
spread out the Ha recombination photons into regions where 
there is no ionizing source. 

5.1. Ionizing and heating sources 



The present scenario is able to explain the line ratios ob- 
served in the WIM much more easily than or equally as well 
as the photoionization models. The higher [S II] A6716/Ha 
and [N II] A6583/Ha line ratios observed in the diffuse Ha 
regions compared to those in H II regions have been be- 
lieved to be the strongest evidence against the theory of a 
dust-scattering origin of the diffuse Ha emission (Reynolds 
1990b). However, we have to note that the argument has been 
based solely on comparison with the line ratios from the bright 
H II regions. If the line ratios from the ionized regions due to 
the late-type OB stars are compared with the line ratios in the 
WIM, the present scenario seems to explain the trend equally 
as well as the previous photoionization models do. 

As already noted in Section 3, mean values of [S II]/Ha and 
[N II]/Ha intensity ratios in the diffuse Ha regions are well 
reproduced with late OB stars (Sembach et al. 2000). Fig- 
ure 3 in Reynolds (1988) shows that the [S II]/Ha ratios ob- 
served in H II regions around late-type stars are higher than 
those found around earlier-type stars. Recently, O'Dell et al. 
(201 1) pointed out that the similarity of the line ratios between 
the Barnard's Loop, the Orion-Eridanus Bubble, and the typ- 
ical WIM samples is striking and probably indicates a com- 
mon set of their physical conditions. The ionizing stars of the 
Barnard's Loop are 09.5 and 09 types. They found that the 
photoionization models with the ionization parameter log U = 
-3.07 ~ -3.67 and the stellar temperature T slsx = 31,000 K - 
40,000 K fully enclose the space of [S II]/Ha versus [N II]/Ha 
occupied by Barnard's Loop, the Orion-Eridaunus Bubble, 
and the WIM samples. It is also noticeable that most of the 
WIM samples in their Figure 7 occupies the region defined by 
a stellar temperature lower than 40,000 K, probably ~ 35,000 
K. Another interesting fact to note is that the elemental abun- 
dances required to explain the line ratios in the WIM are close 
to the B star abundances both in Sembach et al. (2000) and in 
O'Dell et al. (2011). These results strongly support the late 
OB star origin of the diffuse Ha emission. 

Weak [O III] A5007/Ha and He I A5876/Ha emission line 
ratios also indicate that the spectrum of the diffuse interstel- 
lar radiation field is significantly softer than that from the 
average Galactic O star population (Reynolds & Tufte 1995; 
Madsen et al. 2006). 

The increase of the line ratios [S II]/Ha and [N II]/Ha 
can be understood as a significant ioninzation contribution 
by late OB stars to the WIM. However, their highest values 
observed in the WIM may need other non-ionizing heating 
sources, such as shocks, photoelectric heating, and/or turbu- 
lent mixing layers (Reynolds & Cox 1992; Slavin et al. 1993; 
Collins & Rand 2001). Other heating sources may include 
cooling, falling galactic fountain gas initially raised from the 
midplane by supernovae and microflares from magnetic re- 
connection (Raymond 1992; Shapiro & Benjamin 1993). It 
should be noted that the pure photoionization models by O- 
type stars in the midplane also need additional heating sources 
to explain the highest ratios (Reynolds et al. 1999). The emis- 
sion lines due to the heating sources that provide the high- 
est line ratios would also be scattered into more extended re- 
gions than those originally produced by the heating/ionizing 
sources. Therefore, the dust-scattering will hamper the clear 
identification of the original ionization sources. 

As noted by Hoopes & Walterbos (2003), the [S II]/[N II] 
ratio is almost independent of ionizing stellar temperature 
but depends on the ionization parameter. The photoioniza- 
tion models explain the enhanced [S II]/Ha line ratios by a 
diluted radiation field and thus the lowered ionization param- 
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eter. Hence the models cannot reproduce the observed con- 
stancy of the [S II]/[N II] ratios. Instead, the dust-scattering 
scenario can naturally explain the constancy of the line ratios. 

5.2. Physical properties 

Dust scattering has a significant impact on line ratios which 
are used to determine physical properties of H II regions 
(O'Dell & Goss 2009; O'Dell & Harris 2010). Scattered light 
is enhanced at shorter wavelengths, which can lead to ob- 
served Balmer line ratios that are theoretically impossible. 
This effect can also lead to overestimates of the electron tem- 
peratures derived from the auroral and nebular line ratios of 
forbidden lines, such as [N II] A5755 /[N II] A6583 which has 
been used to derive the temperature of the WIM. Serious dis- 
crepancies between the Ha intensity and radio and/or IR in- 
tensities could be found when a large portion of the Ha in- 
tensity is due to dust-scattering. This effect was found in the 
Extended Orion Nebula (EON). O'Dell & Goss (2009) com- 
bined 327.5 MHz radio observations and optical spectroscopy 
to study conditions in the EON. They found an increase in the 
ratio of emission measures derived from the H/3 line and the 
327.5 MHz radio continuum with increasing distance from 
the dominant photoionizing star 8 l Ori C, indicating the in- 
creasing contribution of a dust-scattering component to the 
H/3 intensity with the distance. 

Recent observations using the Wilkinson Microwave 
Anisotropy Probe have found that the ratio of the free-free 
radio continuum to Ha is surprisingly low in the WIM 
(Davies et al. 2006; Dobler & Finkbeiner 2008; Dobler et al. 
2009; Gold et al. 2011). Dong & Draine (2011) proposed a 
three-component model consisting of a mix of (1) hot gas 
currently being photoionized, (2) gas that is recombining and 
cooling after removal of a photoionizing source, and (3) cold 
H I gas. In the standard model for explaining the observed 
intensity ratios of the free-free radio continuum to Ha and 
[N II] A6583 to Ha, they assumed that the scattered fraction 
of the Ha originating from the hot gas is 20%. The frac- 
tions of the hot, cooling, and H I gases were found to be 22, 
56, and 2%, respectively. The model predicted that the pho- 
toionization should switch off and the gas begins to cool and 
recombine in a much shorter time (1 x 10 5 yr) than O star 
lifetimes (~ 3 x 10 6 yr). Such a short photoionization time 
scale is rather surprising in that O stars have been generally 
favored as the source of ionization for the WIM. They there- 
fore suggested that the ionizing radiation for the WIM may be 
provided in large part by runaway O and B stars, with space 
velocities > 100 km s" 1 . In the analysis, they assumed the 
[N II]/Ha line ratio of ~ 0.4 for the WIM, which is well re- 
produced by late OB-type stars (Sembach et al. 2000). The 
ratio of free-free to Ha is estimated to be tp = 0.16 kJy sr -1 
R" 1 for the WIM with a temperature of 8000 K (Equation (11) 
in Dong & Draine 201 1), which is about 1.9 times larger than 
the observed ratio of 0.085 kJy sr" 1 R" 1 . If we assume that 
^37% of the observed Ha intensity is from dust-scattering 
and the temperature of the WIM is 8000 K, the fractions of 
the hot and cooling gases are found to be ~44% and ~48%, 
respectively. 

The WIM is thought to have a temperature 2000 K higher 
than the bright classical H II regions (Reynolds et al. 2001; 
Madsenetal. 2006). However, O'Dell et al. (2011) found 
no evidence that the WIM components have a systematically 
higher temperature than the well known higher-density H II 
regions. As already noted, the temperature derived from 



[N II] A5755 /[N II] A6583 may be uncertain. Tempera- 
ture of the Barnard's Loop is found to be <~ 5960-6100 K 
(Heiles et al. 2000; O'Dell et al. 2011). From the similarity 
of the physical conditions between the Barnard's Loop and 
the typical WIM sample, we may assume an extreme case 
that the temperature of the WIM is ~ 6100 K, and there is 
no contribution from either cooling or cold H I gas. The as- 
sumptions may be rather extreme, but still plausible at some 
sightlines since the physical conditions and constribution of 
dust scattering in individual sightlines would generally dif- 
fer from the average and vary from sightline to sightline 
(Madsen et al. 2006; Witt et al. 2010). The free-free to Ha 
ratio of (1-0.37)^(6100 K) ~ 0.085 is then obtained, which 
agrees with the observed ratio. Therefore, relatively late O 
and/or early B-type stars, as proposed in the present study, 
appear to explain the observed values of not only optical line 
ratios but also free-free to Ha intensity ratios even without 
considering a significant contribution of either cooling or cold 
H I gas. However, we emphasize that the aim of this study is 
to explain averaged observational results and our results do 
not indicate that all the diffuse Ha emission is produced by 
late OB stars outside bright H II regions and dust scattering. 
Large scatter found in Figure 3 also indicates that not only the 
physical conditions and contributions of various components 
(dust-scattered component, hot and cooling gases) but also the 
ionizing and/or heating sources of the WIM vary significantly 
from sightline to sightline. 

Low extinction at 2.17 pm Br7 emission line should allow 
to probe the WIM throughout the Galactic plane. A pilot 
survey of the Br7 emission carried out using the Goddard- 
Wisconsin near-IR cryogenic spectrometer revealed that the 
volume filling fraction of the Br7-emitting medium is typi- 
cally around 1%, indicating the Br7 emission is likely related 
to more compact sources (Kutyrev et al. 2001, 2004). If the 
Lyc photons are indeed leaked out of the bright H II regions in 
the Galactic plane and photoionize the surrounding medium 
to produce the diffuse Ha emission, as usually supposed, the 
filling fraction of the Br7-emitting gas in the Galactic plane 
should be equal to or higher than the filling fraction of ~10% 
obtained by the Ha observations (Kulkarni & Heiles 1987; 
Reynolds 1991; Peterson & Webber 2002). This discrepancy 
strongly suggests that the Lyc leaked out of bright H II regions 
is unlikely the major source of the diffuse Ha background and 
a significant portion of the Ha background originates from 
dust-scattering. More detailed observations using hydrogen 
recombination lines in the near-IR wavelenths, such as Paa 
(1.87 (im), Bra (4.05 /im), and Br7, would be needed to un- 
derstand the origin of the diffuse Ha emission. The Multi- 
purpose Infra-Red Imaging System (MIRIS) is being devel- 
oped to survey the Galactic plane at Paa emission line, which 
may help to understand the origin of the diffuse Ha emission 
(Hanetal. 2010). 

The volume filling fraction and scaleheight of the WIM 
were estimated from the assumption that the pulsar dispersion 
measure and the Ha photons probe the same ionized medium. 
However, as a large portion of the diffuse Ha background at 
high latitudes may be the result of dust-scattering, the filling 
fraction and scaleheight must be reexamined. We also found 
that some of the small-scale features in the FUV sky coincided 
with the Ha sky, but not always. A detailed examination of 
small-scale correlation between the Ha and FUV background, 
which should provide better understanding of the origin of the 
diffuse Ha emission, will be presented in the future. 



Comparison of the Ha and FUV continuum backgrounds 
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5.3. Radiative transfer mechanisms 

Assuming the conventional photoionization models of the 
Ha background, the ISM may have low-density paths and 
voids that allow for ionizing photons from midplane OB stars 
to reach and ionize gas many kiloparsecs above the Galactic 
plane (Wood et al. 2010). Pathways that provide lower than 
average densities to high latitudes could enhance both scat- 
tered FUV and gas ionized by the Lyc from the same OB 
stars in the Galactic plane. Stars that can reproduce most 
of the FUV continuum background in the radiative transfer 
models (e.g., Witt et al. 1997; Schiminovich et al. 2001) are 
located within only a few hundred pc from the Sun. FUV 
photons are predominantly scattered in the forward direction 
(Witt et al. 1997; Schiminovich et al. 2001; Lee et al. 2008), 
implying that dust-scattering occurs mostly in between the 
observer and the source. Preliminary results from our Monte- 
Carlo 3D radiative transfer models for the FUV continuum 
background indicate that most of the FUV continuum back- 
ground originates from stars and dust located within a vol- 
ume with a size of ±1 kpc centered at the Sun (Seon et al., in 
preparation). The Ha photons of the WIM are thought to orig- 
inate within 2-3 kpc distance (e.g., Reynolds 1984), implying 
that Lyc photons would also originate from a volume within 
a few kpc. In this regard, we note that the dust-scattering 
cross-section at 1550A is ~ 1.3 X 10 4 times lower than the 
photoionization cross-section at 912A (Weingartner & Draine 
2001; Verner & Yakovlev 1995), assuming a Galactic dust- 
to-gas ratio. This suggests that FUV photons should origi- 
nate from much (probably ~ 10 4 times) larger volume than 
Ha photons, which is in contrast with the radiative transfer 
models of the FUV continuum background. Meanwhile, the 
dust-scattering cross-section for Ha photons is lower than that 
for FUV (1550 A) by a factor of ~ 2 (Weingartner & Draine 
2001), indicating that dust-scattering of the FUV photons will 
occur at closer distances than Ha photons. Therefore, the 
dust-scattering of Ha photons provide a much easier way to 
explain the diffuse Ha background than the conventional pho- 
toionization models. 

Wood & Reynolds (1999) found that the Ha intensity from 
midplane H II regions that is scattered by dust at high lati- 
tudes is in the range 5-20% of the total intensity. However, 
we should note that in-situ Ha recombination photons at high 
latitudes would be also scattered and this component was not 
included in the estimation. In Figure 3 of Wood & Reynolds 
(1999), the scattered Ha intensity that is originating from 
the WIM through in-situ recombination is compatible with 
the scattered light from bright H II regions (point sources) at 
some high latitudes, indicating higher fraction of dust scatter- 
ing than 5-20%. Moreover, as noted in Witt et al. (2010), the 
highly structured distribution in the ratio of Ha scattering to 
Ha in situ emission was not accounted for in the model of 
Wood & Reynolds (1999). The dust-scattered light of in-situ 
Ha emission (originating mostly from late OB-type stars) at 
high latitudes, and the complex structures of the ISM and ion- 
izing sources may result in a higher fraction of dust-scattering 
than predicted by the simplistic calculation. 

If late OB stars produce about half of the diffuse Ha inten- 
sity, the dust-scattered fraction of the Ha intensity is more or 
less at the same level as that due to in-situ recombination. The 
direct FUV starlight and the diffuse FUV background were 
found to contribute approximately the same amounts to the 
total FUV intensity on average (Seon et al. 2011). The sim- 
ilarity of the dust-scattered fractions between the FUV con- 



tinuum and Ha backgrounds supports again the idea wherein 
their sources and radiative transfer mechanisms are similar. 

6. SUMMARY 

In summary, we found that the Ha/FUV intensity ratio in- 
creases with the Ha intensity and FUV hardness ratio, imply- 
ing that late OB-type stars mostly produce the diffuse Ha at 
high latitudes. We also found that the Ha correlates with the 
FUV background very well. The correlations strongly suggest 
that a large portion of the diffuse Ha photons originate from 
late OB-type stars and dust-scattering. We therefore propose a 
scenario wherein both late OB-type stars and dust-scattering 
play significant roles in the diffuse Ha emission. The Ha 
background at high Galactic latitudes may mostly originate 
from late OB stars existing outside the bright H II regions and 
dust-scattering of the Ha photons. The proposed scenario ap- 
pears to reproduce most of the observed line ratios, such as av- 
erage values of [S II]/Ha, [N II]/Ha, [O III]/Ha and He 1/Ha, 
and free-free to Ha ratios. Their extreme values may need 
additional heating sources, as is also needed in the pure pho- 
toionization models. Our scenario is also in better accord with 
the relative constancy of the [S II]/[N II] intensity ratio, which 
cannot be explained with the standard photoionization mod- 
els. 
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Figure 1. Mollweide projections of the diffuse Hce to FUV intensity ratio maps. Top panel shows the ratio map before dust-extinction correction. Bottom 
panel shows the ratio map after dust-extinction correction. The dust-extinction correction was performed assuming the uniformly mixed gas and dust. Here, 1 
R (rayleigh) = l0 6 /4n photons cm -2 s _I sr~' and 1 CU (continuum unit) = photons cnr 2 s -1 sr~' A -1 . Galactic coordinates centered at (l,b) = (0°,0°) with 
longitude increasing toward the left are shown with latitude and longitude lines on a 30° grid. 
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Figure 2. Histograms of the Ha/FUV ratios for the H II regions and WIM. Extinction was (a) not corrected and (b) corrected assuming that the emitting gas is 
coextensive with dust. A indicates the difference of the mean intensity ratios in logarithmic scale for various thresholds. The ratios calculated for various stellar 
types are also indicated. 
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Figure 3. (a) The Ha/FUV intensity ratio versus Ha intensity, (b) The FUV Hardness ratio (1370-1520 to 1560-1710A) versus the Ha/FUV intensity ratio. 
Black lines represent the average (Ha/FUV intensity or FUV hardness) ratios for given abscissa values after dust-extinction correction. Standard deviations of 
the Ha/FUV ratios from its mean values are also shown. Red dots represent the ratios before the correction of dust-extinction. 
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Figure 4. Correlation of the diffuse Ha background with (a) IR 100 (im emission, and (b) neutral hydrogen column density N(H I). The left panels show 
two-dimensional histogram and the right panels contours of the histograms. Numbers in the left panels are the correlation coefficients estimated in logarithmic 
scale. The contours correspond to 0.7, 0.3, and 0.1 of the maximum values of the histograms. See Figure 21(c) in Seon et al. (201 1) for the comparison of the Ha 
background with the FUV continuum background. The correlation coefficient of the Ha versus FUV intensities is 0.81 from Figure 21(c) in Seon et al. (2011). 
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Figure 5. /sin \b\ versus sin \b\ for the neutral hydrogen column density N(H I), 100 (im emission, Ha, and FUV intensities. Curves were shifted arbitrarily for 
comparison of the shapes. Top curves represent the variation of JV(H I) and the 100 (im emission in red and black colors, respectively. Middle and bottom curves 
show the variation of the Ha (red) and FUV (black) intensities. Middle curves represent the values after dust-extinction correction and bottom curves before the 
correction. Units are denoted in square brackets. 
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Figure 6. Correlation of the diffuse Ha intensity with neutral hydrogen column density (left), 100 /im emission (middle), and the diffuse FUV intensity (right). 
Red and black dots represent the data points before and after the dust-extinction correction, respectively. Numbers outside and inside the parentheses are 
correlation coefficients estimated after and before the dust-extinction correction, respectively. 



